Coulomb drag in quantum circuits 
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We study drag effect in a system of two electrically isolated quantum point contacts (QPC), coupled by 
Coulomb interactions. Drag current exhibits maxima as a function of QPC gate voltages when the latter are 
tuned to the transitions between quantized conductance plateaus. In the linear regime this behavior is due to 
enhanced electron-hole asymmetry near an opening of a new conductance channel. In the non-linear regime the 
drag current is proportional to the shot noise of the driving circuit, suggesting that the Coulomb drag experiments 
may be a convenient way to measure the quantum shot noise. Remarkably, the transition to the non-linear regime 
may occur at driving voltages substantially smaller than the temperature. 

PACS numbers: 73.63.-b, 73.63.-Rt 



Drag effect in bulk 2D systems is well established exper- 
imentally [1-6] and studied theoretically [7-10]. By now it 
is one of the standard ways to access and measure electron- 
electron scattering. Very recently a number of experiments 
were performed to study Coulomb drag in quantum con- 
fined geometries such as quantum wires [11-14], quantum 
dots [15, 16] or quantum point contacts (QPC) [17]. In these 
systems a source-drain voltage V is applied to generate cur- 
rent in the drive circuit while an induced current (or voltage) 
is measured in the drag circuit. Such a drag current is a func- 
tion of the drive voltage V as well as gate voltages, which 
controls transmission of one or both circuits. Figure la shows 
an example of such a setup, where both drive and drag circuits 
are represented by two QPC's. 

It was reported [11-13, 16, 17] that the drag current ex- 
hibits maxima for specific values of the gate voltage, where 
the drive QPC is tuned to an opening of another conductance 
channel. This observation is depicted schematically in Fig. lb. 
It was attributed to the shot noise of the drive QPC [15-18], 
which is known [19, 20] to exhibit a qualitatively similar be- 
havior. The idea is that the drag circuit serves as a detector and 
a rectifier of the quantum shot noise in the drive circuit. Al- 
though plausible and in a certain regime indeed correct, this 
mechanism differs substantially from the one familiar from 
the bulk 2D drag effect. In the latter case drag may be in- 
terpreted [9] as a rectification of nearly equilibrium classical 
thermal fluctuations in the drive circuit. As a result the drag 
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current is a power-law function of the temperature (~ T in 
many cases [21]). Such a rectification is only possible due to 
electron-hole asymmetry in both circuits (otherwise drag cur- 
rents of electrons and holes cancel each other). In the bulk 
systems the asymmetry is due to a small curvature of the par- 
ticles dispersion relation near the Fermi energy. 

Mesoscopic and quantum circuits with the spatial dimen- 
sions less than the temperature length Lj = Vf/T and volt- 
age length Ly = Vf/eV differ from the bulk 2D systems in 
several important ways, (i) The electron-hole symmetry in 
such devices is broken much stronger than in bulk systems. In 
mesoscopic devices this is due to a random configurations of 
impurities [23], while in the QPC's the effect is due to the en- 
ergy dependence of transmission coefficients. Because of the 



FIG. 1: (Color online) a) Two coupled QPCs and surrounding elec- 
tric circuitry. The Coulomb coupling is due to mutual capacitances 
C c . Gate voltage V g control transmission of e.g. drive QPC. b) 
Schematic representation of linear conductance of the drive QPC 
along with the drag current as a function of the gate voltage. 



latter the hole's transmission probability is typically less than 
that of the electrons, (ii) The spatial inversion symmetry may 
be broken by left-right asymmetry of the circuit design. As 
we show below, this makes two polarities of the drive voltage 
V to be essentially non-equivalent, (iii) Because of the above, 
the quantum circuits may be easily driven out of the linear re- 
sponse domain (unlike the bulk systems). Typically a voltage 
needed to drive a quantum circuit into a nonlinear regime is 
parametrically smaller than the temperature. 

In this paper we study drag effect between two QPC's, 
Fig. la. We assume weak interaction between the two cir- 
cuits mediated by mutual capacitances C c , Fig. la. Since the 
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external circuits typically include also dissipative elements, 
the actual interaction is, in general, frequency-dependent [15] 
and determined by a matrix of trans-impedances (see be- 
low). Because of the weak coupling the drag current I D is 
small and therefore the drag circuit is assumed to be close 
to equilibrium. On the other hand, the drive circuit may be 
substantially out-of-equilibrium, due to an applied bias V. 
With these assumption we evaluate the drag current I D in the 
second-order in the inter-circuit interactions in the framework 
of the Keldysh diagrammatic technique (to account for non- 
equilibrium conditions of the drive circuit). Details of the cal- 
culations are reported as a supplementary material. 

We show that at sufficiently small driving voltage V the drag 
current is linear Id °c V. In this regime the mechanism of the 
drag is similar to that in the bulk 2D systems: i.e. rectification 
of near-equilibrium thermal noise. Consequently Id x T 2 at 
small temperatures. The rectification relies on the electron- 
hole asymmetry, which is due to energy dependence of the 
transmission probability in a given channel. The asymmetry 
is the strongest near an opening of a new conductance chan- 
nel. Indeed, in this case thermally excited electrons are much 
more likely to be transmitted than the holes. Hence the behav- 
ior sketched in Fig. lb (though with no relation to the quan- 
tum shot noise). At larger drive voltages I D °c V 2 and the 
effect is indeed due to the detection of the excess shot noise 
of the drive circuit [15-18]. The energy dependence of the 
transmission probability is not required in this regime, and I D 
is proportional to the celebrated Fano factor [19, 20, 30]. Re- 
markably, the crossover between the two regimes takes place 
at eV ~ T 2 /A <sc T, where A is an energy scale of the trans- 
mission probability. 

Quantitatively we found the following expression for the 
drag current: 



Id(V) = j ^Tr[Z{oj)S,{co,V)Z{-co)t 2 {oS)] . 



(1) 



Here and throughout the paper indexes 1,2 refer to the 
drive and drag circuits, correspondingly. The elements 
Zabiw), a,b = R, L of the trans-impedance matrix Z,{(±i) en- 
code inter-circuit coupling. They are defined as Zab(<^>) = 
8<i> a (a))/dIb(cL>), where the corresponding local fluctuating 
currents I a and voltages O a are indicated in Fig. la. 

In Eq. ( 1 ) the drive circuit is characterized by the excess part 
S" b (to, V) = S a t(co, V) - S ab (co, 0) of current-current correla- 
tion matrix S ah (co, V) = J dt e ioj '((6I a (t)6I b (0)+6I b (0)6I a (t))), 
which is known from the theory of quantum shot noise [19, 24, 
25]. In particular 

S ll(u,V) = -^-Yj J de[B RR {e)\t R n {e + )\ 2 \t R n {e-t 

+J B LL (e)[l -r; L (e + )r^(e_)][l -r; L (e_)r^(6 + )] 

+ B LR (e)\r L n (e + )\ 2 \t R (^)\ 2 + fi SL (e)|t*(e + )| 2 |r„ L (e_)| 2 ] , (2) 

with the similar expressions for S ut, S R l and S RR com- 
ponents, see Ref. [24] and supplementary materials for de- 
tails. Here Rq = ^ is quantum resistance, |t^(e)| 2 = 



|t„(e + eVz,(s))| 2 are transmission probabilities of the drive 
QPCi, labeled by the transverse channel index n; |r^ (R) (e)| 2 = 
1 - |t^ (S) (e)| 2 and V L - V R = V. The statistical factors 
are B ah (e) = f a (e + )[l - f b {e_)} + f b (e.)[l - f a {e + )], with 
fu,K)(e) = f(e + sVl{ R )) being the Fermi distributions of the 
two leads and e± = e + cu/2. 

The drag circuit in Eq. (1) is characterized by the rectifi- 
cation coefficient r2(w) = F2(w)f z of the ac voltage fluctua- 
tions applied to the (near equilibrium) drag QPC2, where f z is 
third Pauli matrix acting in the left-right space. Rectification 
is given by 

r2(w) = -^Z Jrfe[/(e_)-/(e + )][|t„(6 + )| 2 -|t„(e_)| 2 ]. (3) 

Characteristics of the QPC2 enter through its energy- 
dependent transmission probabilities |t„(e)| 2 . This expression 
admits a transparent interpretation: potential fluctuations with 
frequency a>, say on the left of the QPC, create electron-hole 
pairs with energies e ± on the branch of right moving particles. 
Consequently the electrons can pass through the QPC with 
the probability |t„(e+)| 2 , while the holes with the probability 
|t„(e_)| 2 . The difference between the two gives the dc current 
flowing across the QPC. Notice that the energy dependence of 
the transmission probabilities in the drag QPC is crucial [26] 
in order to have the asymmetry between electrons and holes, 
and thus non-zero rectification r2(w). 

Focusing on a single partially open channel in a smooth 
QPC, one may think of the potential barrier across it as being 
practically parabolic. In such a case its transmission probabil- 
ity is given by 



|t(e)| 2 = (exp{( e y, - e)/A 2 ] + l)"' , 



(4) 



where A2 is an energy scale associated with the curvature of 
the parabolic barrier in the QPC2 and gate voltage V g moves 
the top of the barrier relative to the Fermi energy. This form 
of transmission was used to explain QPC conductance quan- 
tization [27] and it turns out to be useful in application to the 
Coulomb drag problem. Inserting (4) into (3) and carrying out 
energy integration, one finds 



2eA 2 
Ro 



In 



1 + 



sinh 2 (w/2A 2 ) 
cosh 2 (eV 9 /2A 2 ) 



(5) 



for T <sc A2. In the other limit, T » A2, one should replace 
A 2 -> T in the Eq. (5). Notice that for small frequency to <sc 
A2 one has F2 ~ to 2 , making the the integral in Eq. (1) to be 
convergent in oj — > region. 

Linear drag regime. For small applied voltages V one ex- 
pects the response current Id to be linear in V. Expanding 
Si(u), V) to the linear order in V, one finds that only diago- 
nal components of the current-current correlation matrix con- 
tribute to the linear response and as a result 



S 1 ((o,V) = V — 
ocj 



1 w 
coth — 

2T 



ri(w) + OiV') , 



(6) 
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where Ti(co) is obtained from Eq. (3) by substituting transmis- 
sion probabilities of QPC2, by that of QPQ . Inserting Eq. (6) 
into Eq. (1) one finds 



4 r 



. a+(co) d 
da> — — — 

ocj 



w 1 



coth — 
2T 



ri(w)r 2 (o>), (7) 



where dimensionless interaction kernel a+(tj) is expressed 
through the components of the trans-impedance matrix as 
a ± (fi>) = [(\Zll\ 2 - ZlrZrl) ± (\Zrr\ 2 - ZlrZrl)]/2R 2 q . 
Derived equation (7) has the same general structure as the one 
for the drag current in bulk 2D systems [9, 10]. Being sym- 
metric with respect 1 «-> 2 permutation, it satisfies Onsager 
relation for the linear response coefficient. 

Assuming the load impedance of the drag circuit to be much 
larger than that of the drive one Z\ «Z2« Rq and the mutual 
capacitance of the two circuits to be small C c «: Cr^ s , see 
Fig. la, one finds for the low frequency limit a> <K {ZiC s )~ l 
of the interaction kernels 



l <2 



8Rl C 2 C 2 D 



2C 2 + 2C L C R + 2C\ 

c 2 -c 2 

L R 



(8) 



For Z\ — > the drive QPC is shorted and the drag circuit is 
insensitive to the fluctuations. Substituting now Eq. (5) into 
Eq. (7), one finds for e.g. low-temperature regime T <sc Ai >2 



In = — 



V a + (0)n 2 T 2 



1 



R Q 6 AiA 2 cosh 2 (eVy2Ai) ' 



(9) 



where we assumed that the gate voltage of QPC 2 is tuned to 
adjust the top of its barrier with the Fermi energy and wrote 
I D as a function of the gate voltage in QPCi. We have also 
assumed that T <K {Z\C s y l to substitute a+(co) by its dc limit 
Eq. (8). The resulting expression exhibits a peak at V g = 
similar to that depicted in Fig. lb. Yet it has nothing to do 
with the shot noise, but rather reflects rectification of near- 
equilibrium thermal fluctuations (hence the factor T 2 ) along 
with the electron-hole asymmetry (hence non-monotonous de- 
pendence on Vg). For monotonously increasing functions 
|t(e)| 2 in both circuits the linear drag is positive (i.e. currents 
flow in the same direction). 

Nonlinear regime. At larger drive voltages drag current 
ceases to be linear in V. Furthermore, contrary to the linear 
response case, § i(oj, V) does not require energy dependence 
of the transmission probabilities and could be evaluated for 
energy independent |t„| 2 (this is a fare assumption for T, eV «: 
Ai). Assuming in addition T «: eV, one finds a celebrated 
expression for the quantum shot noise [19, 25] 



S 1 (a>,V) = 2 



\eV + cu\ + \eV-(o\ 

To 



^|t„| 2 [i-iy 2 ]f . (io) 



Inserting Eq. (10) into Eq. (1), after frequency integration 
bounded by the voltage, one finds for the drag current [28] 



In = 



eV 2 

~&2Rq 



mo) 2iu 2 [i-it»i 2 ]. (id 



Here again we assumed that the detector QPCi is tuned to 
the transition between the plateaus. We also assumed eV <sc 
(ZiCj)' 1 to substitute a-(a>) by its dc value, Eq. (8). One 
should notice that while a+ > 0, the sign of a_ is arbitrary. For 
a completely symmetric circuit a_ = 0, while for extremely 
asymmetric one |ar_| < a+/2. Although we presented deriva- 
tion of Eq. (11) for T <K eV, one may show that it remains 
valid at any temperature as long as T «: min{Ai, (Z\C s y l ). 

Equation (11) indeed shows that the drag current is due to 
the rectification of the quantum shot noise and hence propor- 
tional to the Fano factor [19]. It again exhibits a generic be- 
havior depicted in Fig. lb, but the reason is rather different 
from the similar behavior in the linear regime. The direction 
of the nonlinear drag current is determined by the inversion 
asymmetry of the circuit (through the sign of or_) rather than 
the direction of the drive current. As a result, for a certain 
polarity of the drive voltage, the drag current appears to be 
negative. 

We discuss now a crossover between the two regimes. As- 
suming that for a generic circuit a+ ~ a_ and comparing 
Eqs. (9) and (11) one concludes that the transition from the 
linear to the nonlinear regime takes place at V ~ V* with 



eV = r 2 /A, « T , 



(12) 



for r « Aj. In the opposite limit, T > Ai, the crossover 
voltage is given by the temperature eV* = T. However, for a 
circuit with an almost perfect inversion symmetry, i.e. a_ «c 
a+, the nonlinear regime may be pushed to substantially larger 
voltages. Such a symmetric circuitry is not well suited for 
detection of the quantum shot noise. 

Mesoscopic circuits. One or both circuits may be repre- 
sented by a multichannel quasi- ID (or 2D) mesoscopic sam- 
ple. In this case |t„(e)| 2 = g(e) is a dimensionless (in units 
of Rq 1 ) conductance of the sample as a function of its Fermi 
energy. Such a conductance exhibits universal conductance 
fluctuations (UCF) [29], that is g(e) = g + 5g{e), where g » 1 
is an average conductance and 5g(e) ~ 1 is a sample and 
energy-dependent fluctuating part. The characteristic scale of 
the energy dependence of the fluctuating part is the Thouless 
energy Ejh = TiD/I?, where D is electronic diffusion constant 
and L is the sample size. Employing Eq. (3), one finds that the 
rectification coefficient of a given mesoscopic sample may be 
estimated as 

2 

T(co) ~±J-^-, {T, a>) « E Th . (13) 
Kq h T h 

On the other hand, the nonequilibrium part of the noise corre- 
lator Eq. (10) exhibits a well-defined average value 

S 1 (oj,V) = 2(\eV + co\ + \eV-cj\)^-, (14) 

the coefficient 1/3 is specific to a quasi- ID geometry [30]. 

In the Coulomb drag setup, where both circuits are repre- 
sented by mesoscopic elements, employing Eqs. (1), (7) along 
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with (13), (14), one finds for the drag current (both linear and 
nonlinear) 



E 2 



eV ^ 

a-T—g 

V-Th 



(15) 



where T < E T h- If the load impedance of the drive cir- 
cuit is Zi ~ g~ l , then linear in V term of Eq. (15) is in 
agreement with the corresponding result of Ref. [23]. The 
crossover between linear and nonlinear regimes takes place at 
eV* = a + T 2 /(a-gE T h) which may be much less than both T 
and Ejh- As a result, one may expect drag current to be sub- 
stantially bigger than the linear response prediction already at 
the very modest bias voltage. 

In summary, we have studied Coulomb drag effect in the 
system of two coupled quantum circuits. In the linear regime 
gate voltage induced oscillations of the drag conductance 
originate from the particle-hole asymmetry, which is encoded 
in the energy dependent transmission probabilities of the 
QPC. The drag conductance follows quadratic temperature 
dependence at low temperatures and is peaked at gate volt- 
ages, which correspond to the transition between QPC con- 
ductance plateaus. Beyond the linear regime the magnitude 
of the drag current is proportional to the current shot noise 
generated in the drive QPC. 
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version of Eq. (2). This work was supported by NSF grants 
DMR-0405212 and DMR-0804266. 
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Supplementary materials 

The purpose of this section is to provide technical details 
needed to derive Eq. (1) of the main paper. To this end, we 
describe each point contact of the quantum circuit Fig. la 
as quasi- ID adiabatic constriction connected to two reser- 
voirs (terminals, probes), to be referred to as left (L) and 
right (R). The distribution functions of electrons in the reser- 
voirs of a driven circuit, are Fermi distributions fuR){e) = 
[exp[(e + eVuR))/T] + l] -1 , with source-drain voltage being 
Vl - Vr = V. In the dragged circuit distributions are assumed 
to be equilibrium Fermi functions. Within each QPC elec- 
tron motion is separable into transverse and longitudinal com- 
ponents. Due to the confinement transverse motion is quan- 
tized and we assign quantum number n to label transverse con- 
duction channels with <p„(r ± ) being corresponding transversal 
wave function. The longitudinal motion is describe in terms 
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of the extended scattering states - normalized electron plane 
waves incident from the left 



(k)e- 

ikx 



X 

x ■ 



— oo 
+00 



and right 



R( , , _ <f>n(r±) j e ikx + r n (k)e ikx x -> +00 
U " (K ' r) ~ ^ \ t n {k)e- ikx x^-00 



(16) 



(17) 



onto mesoscopic scattering region. Here k and v are the elec- 
tron wave vector and velocity, t„(k) and r„(k) are channel spe- 
cific transmission and reflection amplitudes. Second quan- 
tized electron field operator is introduced in the standard way 

*(r, t) = J] r) + t)u R n {k, rj\ , (18) 

ilk 

where i(/„ '(k, t) are fermion destruction operators in the left 
and right reservoirs correspondingly. For the future use we 
define also current operator 



/ fl (f) = X M «»'^ fl(fe ' f) ^' (fe '' f) ' 

nk 



(19) 



which has matrix elements M a m , = 
(e/2im)fdr ± [u* n a (k,r)d x u a n ,(k',r) - d x [K°(k, r)]u a n , (If, r)], 
constructed from the scattering states (16)— (17). Based on 
the orthogonality condition of transverse wave functions 
J dr ± <f> n (r ± )(f>*,(r ± ) = 5„„>, direct calculation gives 



M 



-e6„ 



K(k)r n (k') - 1 r* n (k)t n (k') 
t*„(k)r n (k') t* n (k)t„(k') 



(20) 



and a similar result for M R n , . In Eq. (20) we have suppressed 
phase factors e ±l ( k - k "> x « l, since \k - k'\ ~ L r ' <K x ', and 
the coordinate x is confined by the sample size L «: L7-. On 
the other hand, we do not keep fast oscillating factors e ±2lkFX , 
since x ~ L » k F l . However, one must keep explicitly mo- 
mentum (or equivalently energy) dependence for transmission 
amplitudes, which translates later into particle-hole asymme- 
try factor T(aj). 

Dynamics of if/„ '(k, t) operators is governed by the action 

iSo = f £ ^ W-Vjnik, t) (21) 

^ c ja nk 

defined along the Keldysh contour, where G _1 = (id t - %h + 
eV a ) is Green's function operator with being electron en- 
ergy. Additional subscript j in Eq. (21) labels drive (J = 1) 
and drag ( j = 2) QPC's. As usual for Keldysh technique one 
splits time integration into forward and backward pathes and 
replicates each fermion field into two components if/ => if/* 1 , 
which belong now to the different contour branches. It is con- 
venient also to perform Keldysh rotation 



L 



1 



-1 
1 

(22) 
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FIG. 2: Drag current I D in the second order in inter-circuit inter- 
actions K = Z/(*> (wavy lines). The drag circuit is represented by 
triangular rectification vertex T 2 {oS), while the drive circuit by the 
non-equilibrium current-current correlator Si(a>, V) (loop). 



where & z is third Pauli matrix acting in the Keldysh space. In 
this rotated basis quadratic action (21) gives following elec- 
tron correlators 



-i«<C(*. # % ( k > f ')» = Sij8abGZ(t, f) , (23) 

where Keldysh Green's function matrix has familiar triangular 
structure 



a/3. 



Gt a (t,t')j ■ 



(24) 



Retarded/Advanced/Keldysh components of G"^ are given by 



iG R Ui {t,l') = e(t-t')e- i( ? M(t - f \ 
iGl(t, t') = -6(t' - t)e- i&+eV " )( <- n , 
iGf a (tj') = [1 -2f a (e k )]e-^ +eVM -''K 



(25) 



Having described quantum point contacts individually we 
introduce now the interaction term between them 



sir 

aba/3 JJ - a 



dtdt' Il(t)Kf h (t ■ 



(26) 



Here IjR(L)(f) are current operators (19), on the right (left) 
of the QPCy, coupled by the kernel K(f - f')> which en- 
codes electromagnetic environment of the circuit. Interac- 
tion kernel retarded and advanced components are directly re- 
lated to the trans-impedance matrix of the circuit K ^ \io) = 
Z R< j^\to)/(oJ ± i0), while Keldysh component can be restored 
from the fluctuation-dissipation theorem: K^(o>) = [K R b (a))- 
K A , (a/)] coth(<jj/2r), i.e. we assume the surrounding electric 
environment to be close to equilibrium. 
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Within this formalism, drag current is found by averaging 
I2 over the fermionic degrees of freedom 

I D = jD[^]Tr[^ 2 M^ 2 ] exp(»5[^r]), (27) 

where S = So + Sim and we used expression (19) for the 
current operator. Being interested in the leading perturbative 
result, we expand exp(;»S) to the second order in interaction 
term <Sj nt . This way one obtains 

!d = \§ D[#] Tr[0- 2 M<A 2 ] 

Tr [IJU2] Tr [hKI 2 ] exp (iS Q [^]). (28) 

Remaining Gaussian integral may be evaluating using the 
Wick's theorem. One inserts expressions (19) for the current 
operators into the traces of (28) and takes into the account all 
the possible Wick's contraction between i/r-fields. The latter 
are given by the Green's functions (23). This way we find our 
main result for the drag current [Eq. (1)] shown diagrammati- 
cally in the Fig. 2. The interaction kernels a ± {cS), rectification 
coefficient T 2 (w) and noise power S a t(d>, V) are given explic- 
itly by the following Keldysh traces: 

= -^Tr [Z R (oj)t Q(z) Z A (co)t z ] , (29) 

Q 

S ab ((o, V) = Tr [Gfr x M a Gcr x M h ] , (30) 
T 2 (w) = Tr [G^MGa-oMG^oM] , (31) 

where t and & are two sets of Pauli matrices in left-right and 
Keldysh spaces, correspondingly. 

Remaining steps of algebra concern calculation of the 
traces in Eqs. (29)-(31). For the interaction coefficients a± 
it is straightforward. Let us demonstrate how Eq. (2) for S ll 
is recovered from Eq. (30). To this end we insert current ma- 
trix M L [Eq. (20)] along with G = diag{G L , G R ] into Eq. (30) 
and calculate trace over left-right subspace, which gives 



S LL (a>, V) = e 



An 2 



(32) 



+Tr[G L {e + ,t)fr x G R {e-,Ov x }\r n {t)\ 2 \tn{n 2 
+Tr[G«(e + ^)<T A .G L (e_,r)^]|t„©| 2 |r n (r)l 2 
+Tr[G R ( e+ ,^G R (e^O^]\t„(0\ 2 \tn(O\ 2 ] ■ 

Recall here that G L ( R) are still matrices in the Keldysh sub- 
space. Using Eq. (24) one calculates remaining traces over 
the Keldysh subspace 

Tr [G a fr x G b fr x \ = G R a G A h + G A a G R h + G K a G K h , (33) 



and performs final integration with the help of Green's func- 
tions G R(A \e,& = (e - £ + eV a ± jO)" 1 and G*(e,£) = 
-2nid(e - £, + eV a )[l - 2f a (g)], which follows from Fourier 
transforms of Eq. (25). It is not difficult to see now that each 
Keldysh trace in Eq. (32) defines statistical occupation factors 
B ah (e) used in Eq. (2), namely Tr [GaO-xG&o-^] oc B ab (e). As 
a result, collecting all the factors, one finds from Eq. (33) the 
final form of noise power given by Eq. (2). In complete anal- 
ogy one may calculate S RR and S LR components of the noise 
power: 

Srr(oj,V) = ^-Yj J de[B LL {e)\t L n {e + )\ 2 \t L n {eJ)\ 2 (34) 

+B RR {e)[\ - rf (6 + )r«(e_)][rf(6_)rf(e + ) - 1] 
+J B LR (e)|t^(6 + )| 2 |r^-)| 2 + B fi L(e)|rf (e + )| 2 |t^(e_)| 2 ] , 



(35) 



[B LL ( £ )[t: L ( £+ )t L „(e-) - t; L (6 + )t^(6_)r; L (e_)r^(e + )] 

+B RR (e)[t R (e + )tf (e_) - tf (e + )t«(6_)r«(e + )rf (e_)] 
+B LR (eX\e + )t R (e-)rf{eS)r L n (e + ) 

+B R deK R (e + )t L n (e-)r:He-)r R (e+)] ■ 

Notice that cross-correlation component is negative. 

Finding F 2 (w) one uses the fact that Green's function is di- 
agonal unity in the left-right subspace G = Gfo and faces 
Keldysh trace of the kind 

Tr [g& x G& G& ] = ^ [G R (e)G R (e ± u)G K (e) 
± 

+ G R (e)G K (e ± a>)G A (e) + G K (e)G A (e ± co)G A (e)] . (36) 

To simplify this equation one should decompose each 
Keldysh component of the Green's function using fluctuation- 
dissipation relation G K {e) = [G R (e) - G A (e)][l - 2f(e)] and 
keep in the resulting expression only those terms, which have 
different causality. Combinations having three Green's func- 
tions of the same kind, like G A G A G A and G R G R G R , will not 
contribute. This way, one finds that Tr ^cr J Gd'oGo"o] <* 
[f(e_) - f(e+)], with e± = e ± (o/2. Remaining trace over the 
current vertex matrices M 2 reduces to the transmission proba- 
bilities at shifted energies, namely Tr[M 2 M 2 M 2 ] oc +|t„(e ± )| 2 . 
As the result, imposing remaining e integration and summa- 
tion over transverse channels, one arrives at T 2 (co) in the form 
of Eq. (4) of the main text. This completes our derivation. 



